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The  effective  binary  diffusivity  of  H2  and  H20  in  a  Ni  and  yittria-stabilized  zirconia  (YSZ)  anode  of  the  solid 
oxide  fuel  cells  (SOFCs)  was  measured  between  650  and  800  °C  using  an  electrochemical  cell  consisting 
of  an  oxygen  pump,  an  oxygen  sensor,  and  a  porous  SOFC  anode  pellet.  The  effective  binary  diffusivity 
was  obtained  from  the  relationship  between  the  current  density  across  the  oxygen  pump,  and  the  H2 
partial  pressure  gradient  across  the  anode  sample  measured  using  the  oxygen  sensor.  The  anode  limiting 
current  density  and  concentration  polarization  were  estimated  using  the  experimental  results. 

©  2009  Elsevier  B.V.  All  rights  reserved. 


1.  Introduction 

Solid  oxide  fuel  cells  (SOFCs)  have  gained  significant  interest  due 
to  their  high  energy  conversion  efficiency,  low  pollution  emission, 
modularity,  and  fuel  flexibility.  Among  the  basic  designs  of  SOFCs, 
the  anode-supported  SOFCs  have  been  extensively  investigated 
because  they  exhibit  higher  cell  performance  compared  to  either 
the  electrolyte-supported  or  cathode-supported  designs  [1,2].  In 
the  anode-supported  SOFCs,  the  anode  concentration  polarization 
is  one  of  the  significant  contributions  to  the  performance  loss, 
because  the  thick  anode  impedes  the  rapid  transport  of  the  gaseous 
reactant  and  product  species  [3,4].  In  a  thick  and  porous  anode, 
the  gas  transport  is  mainly  through  a  diffusion  process  [5],  and 
the  effective  diffusivity  of  gas  species  is  determined  by  the  anode 
parameters  such  as  the  pore  size,  volume  fraction  of  pores,  and 
tortuosity.  Therefore,  it  is  important  to  measure  the  effective  diffu¬ 
sivity  of  gases  under  the  operating  conditions  of  SOFCs  in  order 
to  effectively  engineer  the  anode  microstructure  and  lower  the 
concentration  polarization.  Between  30  and  300  °C,  the  effective 
diffusivity  of  mixed  gases  has  been  measured  using  techniques 
such  as  gas  chromatography  [6,7],  thermogravimetry  [8],  and  zero 
length  column  method  [9].  Although  the  effective  diffusivity  of 
02  and  N2  in  the  SOFC  cathode  has  been  measured  using  an 
electrochemical  cell  between  650  and  800  °C  [10],  experimental 
measurements  of  the  effective  binary  diffusivity  of  H2  and  H20  in 
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the  SOFC  anode  have  not  been  reported  to  the  authors’  knowledge. 
In  the  present  work,  the  effective  binary  diffusivity  of  FI2  and  FI2  O  in 
a  porous  Ni  and  yittria-stabilized  zirconia  (YSZ)  cermet  anode  was 
measured  using  an  electrochemical  cell  comprised  of  an  oxygen 
pump,  an  oxygen  sensor,  and  a  porous  anode  sample.  The  diffusivity 
values  between  650  and  800  °C  were  obtained  from  the  relationship 
between  the  current  density  across  the  oxygen  pump  and  the  FI2 
partial  pressure  gradient  across  the  porous  sample  measured  using 
the  oxygen  sensor.  The  temperature  dependence  of  the  effective 
diffusivity  has  been  discussed,  and  the  anode  limiting  current  den¬ 
sity  and  concentration  polarization  have  been  estimated  using  the 
experimental  results. 

2.  Experimental 

Fig.  la  shows  the  schematic  of  the  electrochemical  cell  to  mea¬ 
sure  the  effective  binary  diffusivity  of  FI2-FI20  through  a  porous 
Ni-YSZ  anode.  It  consists  of  a  dense  6mol%  YSZ  tube  3.9  in.  long, 
0.7  in.  in  diameter,  with  a  wall  thickness  of  0.1  in.  obtained  from 
McDaniel  Advanced  Ceramic  Technologies,  with  a  dense  YSZ  disc 
on  one  end  and  a  porous  Ni-YSZ  disc  on  the  other  end.  A  dense 
YSZ  disc  was  prepared  by  uniaxial  die-pressing  of  8  mol%  YSZ  pow¬ 
der  (Tosoh  Corp.)  followed  by  sintering  at  1400  °C  for  8  h.  A  paste 
of  NiO-YSZ  was  prepared  by  mixing  NiO  (J.T.  Baker)  and  8  mol% 
YSZ  powder  with  alpha-terpineol  and  binder  in  a  desired  ratio,  and 
painted  on  both  sides  of  the  YSZ  disc  and  tube  to  form  the  electrodes 
for  the  oxygen  pump  and  sensor,  respectively.  A  porous  Ni-YSZ  disc 
(50vol%  Ni  +  50vol%  YSZ)  was  fabricated  by  the  high  shear  com¬ 
paction  process  [11,12]  as  described  in  our  earlier  work  using  NiO 
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Fig.  1.  (a)  A  schematic  of  the  electrochemical  cell  to  measure  the  effective  binary 
diffusivity  of  H2  and  H20  in  porous  Ni-YSZ  sample,  (b)  A  schematic  of  Ni-YSZ  sample, 
along  with  partial  pressures  of  hydrogen. 


and  8  mol%  YSZ  powder  [13].  YSZ  paste  prepared  by  mixing  8  mol% 
YSZ  powder  with  a  sintering  aid,  binder,  and  alpha-terpineol  in  a 
desired  ratio  was  painted  on  the  edge  of  the  NiO-YSZ  disc  to  prevent 
leaks  at  the  edges  of  the  sample  and  allow  gas  diffusion  through  a 
well-defined  area.  The  open  area  of  the  porous  sample  available  for 
diffusion  was  1  cm2.  After  drying  in  an  oven,  the  YSZ  disc,  YSZ  tube 
and  NiO-YSZ  disc  were  sintered  at  1300  °C  for  2  h  in  air.  Then,  Ni 
meshes  connected  to  Ni  wires  were  attached  to  the  electrodes  on 
the  YSZ  disc  and  the  YSZ  tube  using  a  NiO-based  adhesive  prepared 
with  NiO  powders  and  Pyro  Putty  2400-T  (Aremco  Products  Inc.). 
The  YSZ  disc  and  Ni-YSZ  disc  were  attached  to  the  two  ends  of  the 
YSZ  tube  using  a  glass  paste  (FuelCellStore).  The  entire  setup  was 
fired  at  850  °C  under  a  flow  of  300  cm3  min-1  of  humidified  form¬ 
ing  gas  (5%  H2  +  95%  Ar)  containing  3%  H20  for  4h  to  convert  NiO 
to  Ni  and  obtain  a  leak-tight  seal. 

The  electrochemical  cell  was  place  in  the  furnace,  and 
measurements  were  performed  between  650  and  800  °C  with 
150  cm3  min-1  of  humidified  H2  containing  3%  H20  using  a  model 
263  A  potentiostat  from  Princeton  Applied  Research.  A  constant 
current  in  the  range  0.01-0.05  A  was  applied  through  the  YSZ  disc 
to  pump  the  oxygen  into  the  tube,  and  the  Nernst  potential  of  the 
oxygen  sensor  was  continuously  measured  until  a  steady  state  was 
established.  A  steady  state  is  assumed  to  have  been  reached  when 
the  EMF  measured  by  the  oxygen  sensor  does  not  change  apprecia¬ 
bly  with  time. 

After  testing,  the  porous  Ni-YSZ  sample  was  sectioned  and 
impregnated  with  epoxy  in  vacuum.  After  the  epoxy  hardened,  it 
was  polished  down  to  a  1  p,m  finish,  and  the  microstructure  and 
porosity  of  its  cross-section  was  examined  using  scanning  electron 
microscopy  (SEM). 


3.  Theoretical  analysis 

As  oxygen  is  pumped  into  the  YSZ  tube,  the  partial  pressure  of 
H2  inside  the  tube  decreases,  and  that  of  H20  increases.  The  par¬ 
tial  pressure  gradients  across  the  porous  Ni-YSZ  sample  generates 
fluxes  of  H2  C/h2)  and  H20  C/h2o)  into  and  out  of  the  tube,  respec¬ 
tively,  and  in  a  steady  state,  the  fluxes  of  H2  and  H20  across  the 


porous  Ni-YSZ  sample  are  related  to  the  flux  of  02  (Jo2 )  across  the 
YSZ  disc  by  the  current  density  in  the  oxygen  pump  (i)  as: 


Jh2  =Jh2o  =  2/q2  =  2^ 


(1) 


where  Fis  the  Faraday  constant.  Since  the  total  pressure  is  fixed  and 
the  fluxes  of  H2  and  H20  are  in  opposite  directions  C/h2  +  Jh2o  =  0), 
Jh2  and  Jh2o  can  be  expressed  assuming  one-dimensional  partial 
pressure  gradients  as: 


Jh2  2p  \-H2oVnH2 


dh2-h2o  dpn2 


dx 


(2) 


Jh20  -  nc  -  ^4j-H20^nH20  - 
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RT 
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where  D^_H20  is  the  effective  binary  diffusivity  of  H2  and  H20  in 
the  porous  Ni-YSZ  sample,  nH2  and  nH2o  are  the  concentrations  of 
H2  and  H20,  respectively,  R  is  the  gas  constant,  T  is  the  temper¬ 
ature,  and  Ph2  and  pH2o  are  the  partial  pressures  of  H2  and  H20, 
respectively.  Rearranging  and  integrating  Eqs.  (2)  and  (3)  across 
the  porous  Ni-YSZ  sample: 


[p*2  fl  RTi 
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(5) 


where  and  p°  are  the  H2  partial  pressures  inside  and  outside 

the  tube,  respectively,  p^0  and  p^Q  are  the  H20  partial  pressures 
inside  and  outside  the  tube,  respectively,  and  l  is  the  thickness  of  the 
porous  Ni-YSZ  sample.  Fig.  lb  shows  a  schematic  of  a  Ni-YSZ  sam¬ 
ple,  with  a  thickness  and  the  partial  pressures  of  hydrogen  labeled. 
From  Eqs.  (4)  and  (5)  and  Fig.  lb,  the  partial  pressures  of  H2  and 
H20  inside  the  tube  at  a  fixed  current  density  (galvanostatic  mode) 
can  be  expressed  as: 


„(*) 


Ph2=P  h2- 


RTl 


2FDZ-h20 


m  0  RTI 

Ph2o=Ph2o  +  : - - - 1 


2FD 


,  eff 
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(7) 


In  this  experiment,  p°2  and  p^Q  are  fixed  (97%  H2  +  3%  H20), 

and  p[^  and  p^0  can  be  determined  from  the  measured  Nernst 
potential  using  the  oxygen  sensor  on  YSZ  tube.  The  Nernst  potential, 
F,  is  given  by 


(8) 


where  pQ  and  p°  2  are  the  02  partial  pressures  inside  and  outside 
the  tube,  respectively.  Thus,  the  sensor  EMF  provides  a  direct  mea¬ 
sure  of  the  concentration  polarization  across  the  anode  at  a  fixed 
current  density.  Once  the  p[£  is  calculated  at  a  fixed  current  den¬ 
sity  from  the  sensor  EMF  using  Eq.  (8),  p^  andp^0  can  be  obtained 
using  the  equilibrium  constant  (I<eq)  for  the  equation  2H2  (g)  +  02 
(g)  =  2H20  (g)  (given  below  in  Eq.  (9),  and  the  fact  that  ptotai  =  1  atm 
^Ph2  +  Ph2o  everywhere  within  the  anode  and  inside  the  YSZ  tube; 
specifically  inside  the  YSZ  tube,  p^0  =  1  -  p(^. 


I<eq  = 


ph2o 

P»2p°2 


(9) 
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Fig.  2.  SEM  image  of  Ni-YSZ  anode  sample. 


Time  (sec) 

Fig.  3.  The  Nernst  Potential  across  the  oxygen  sensor  as  a  function  of  time  at  800  °C 
with  the  applied  current  density  of  0.05  A  cm-2. 


Therefore,  £^_h2o  can  be  obtained  from  the  slope  of  either  the 
p^  vs.  i  or  the  p[^0  vs.  i  plot  using  either  Eq.  (6)  or  (7). 


4.  Results  and  discussion 

Fig.  2  shows  the  microstructure  of  porous  Ni-YSZ  sample.  The 
thickness  was  measured  to  be  750  p,m,  and  the  porosity  to  be  22%. 
The  mean  pore  size  was  measured  to  be  1.3  p,m. 

Fig.  3  shows  typical  plot  of  the  measured  Nernst  potentials  as  a 
function  of  time  at  800  °C  at  a  fixed  current  density  of  0.05  A  cm-2. 
Similar  measurements  were  made  at  800,  750,  700,  and  650  °C 
between  0.01  and  0.05  A  cm-2 .  The  steady  state  p[£  ’s  were  obtained 
using  the  measured  Nernst  potentials  and  Eqs.  (8)  and  (9)  and  Fig.  4 
shows  the  p[^  vs.  current  density  plot  at  800  °C.  The  slopes  at  differ¬ 
ent  temperatures  were  obtained  by  fitting  into  the  linear  equation, 
and  £>h{-h2o’s  were  calculated  using  Eq.  (6).  Dj^_H  0  values  were 
0.070  cm2  s-1  at  800  °C,  0.060  cm2  s-1  at  750  °C,  0.055  cm2  s”1  at 
700  °C,  and  0.046  cm2  s_1  at  650  °C.  In  previous  studies,  De^_H^Q 
in  the  Ni-YSZ  anode  obtained  by  fitting  i-V  curves  to  the  polar¬ 
ization  model  has  been  reported  to  be  0.07-0.82  cm2  s-1  at  800  °C 
in  anodes  with  volume  percent  porosity  in  the  range  of  23-76% 
[5,13-16].  In  the  phenomenological  theory  of  diffusion  through 
porous  bodies,  the  A-B  effective  binary  diffusivity,  De£_B,  is  related 
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Fig.  4.  The  H2  partial  pressure  inside  the  tube  (p[J  )  as  a  function  of  current  density 
at  800  °C. 


to  the  A-B  binary  diffusivity,  DA_B,  by 

Df-B  =  TDa“B  (10) 

where  Vv  is  the  volume  percent  of  porosity,  and  r  is  the  tortuos¬ 
ity  factor  [17].  Eq.  (10)  shows  that  the  effective  binary  diffusivity 
is  a  function  of  porosity  and  tortuosity  in  the  anode,  which  was 
systematically  studied  in  our  previous  study  [14].  The  diffusivities 
measured  here  are  smaller  than  those  we  have  reported  previ¬ 
ously,  presumably  because  of  the  lower  volume  percent  porosity 
in  the  present  anode  sample.  The  Chapman-Enskogg  model  gives 
the  Da_b  as  [17]: 


f^A-B  = 


0.001 86  r3/2((l/MA)  +  (1/Mb))1/2 
PfloA-B 


(11) 


where  MA  and  MB  are  the  molecular  weights  of  the  components  A 
and  B,  respectively,  p  is  the  total  pressure,  Q  is  the  collision  inte¬ 
gral,  and  crA_B  is  the  collision  diameter.  Using  the  Q  and  crA_B  data 
from  Cussler  [17],  Dh2-h2o  is  calculated  to  be  7.6  cm2  s-1  at  800  °C 
[17],  and  using  the  measured  volume  porosity  of  the  anode  and  Eq. 
(10),  the  tortuosity  factor  is  estimated  to  be  ~22.  The  tortuosity 
factor  is  typically  between  two  and  six  in  porous  bodies  with  a  rel¬ 
atively  coarse  structure,  and  the  value  of  21  has  been  reported  on 
the  SOFC  anodes  with  a  few  microns  of  pore  size  and  23%  of  volume 
porosity  [5,1 5,1 7-20].  Although  it  is  difficult  to  estimate  it  directly, 
since  the  average  pore  size  of  the  Ni-YSZ  sample  used  in  this  work 
was  1 .3  |jim,  the  value  of  ~22  obtained  in  this  work  is  quite  reason¬ 
able.  Eqs.  (10)  and  (11)  show  that  £>h{-h2o  proportional  to  T3!2, 
and  D$  H  Q  vs.  T3!2  plot  in  Fig.  5  shows  the  linear  dependence  as 
expected. 

The  limiting  current  density  and  concentration  polarization  of 
the  anode  can  be  predicted  using  the  measured  £>h{-h2o  values  to 
estimate  the  contribution  of  the  anode  concentration  polarization 
to  the  total  performance  loss  of  SOFCs.  In  operating  cells,  the  anode 
concentration  polarization  occurs  due  to  the  slow  mass  transport 
of  gas-phase  reactant  (H2)  and  product  (H20)  species  through  the 
porous  Ni-YSZ  anode,  and  can  be  expressed  as  [5,13-15]: 


P  a, cone  — 


RT  M 
— —  ln(  1 
2  F  v 


i  ,  RT , 
ias)+  2  F  n 


JkL) 

Ph2oW 


(12) 


where  ias  is  the  anode  limiting  current  density,  and  p°2  and  pf^Q 
are  the  partial  pressures  of  H2  and  H20  in  the  anode  chamber 
gas  stream,  respectively.  The  anode  limiting  current  density  is 
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the  diffusivity  measured  in  the  present  work  is  smaller  than  those 
of  our  own  estimates  from  curve  fitting  voltage-current  density 
curves  of  single  SOFCs  made  in  our  lab  [14].  The  anode  concentra¬ 
tion  polarization  can  be  reduced  by  increasing  the  anode  porosity 
and  decreasing  the  thickness.  The  methodology  presented  in  this 
paper  can  be  used  to  characterize  the  gas  transport  through  the 
anode  and  engineer  the  microstructure.  We  plan  to  conduct  studies 
on  other  SOFC  anode  materials  using  this  technique  and  to  optimize 
the  anode  microstructure  and  porosity. 

5.  Conclusion 
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Fig.  5.  A  plot  of  £>h2-h2o  vs-  73/2  between  650  and  800  °C. 


Fig.  6.  Anode  concentration  polarization  as  a  function  of  current  density  at  700  and 
800 °C  calculated  using  Eq.  (12). 


the  current  density  at  which  the  hydrogen  partial  pressure  at  the 
anode-electrolyte  interface  becomes  zero,  and  given  by  [5,13-15]: 


ias 


2  FP°h2D 


eff 

H2-H20 


RTla 


(13) 


where  la  is  the  thickness  of  the  anode.  The  anode  limiting  current 
density  can  be  predicted  from  experimentally  measured  £>h{-h2o 
values  using  Eq.  (13),  and  they  are  calculated  to  be  1.98  A  cm-2  at 
800  °C,  1 .78  A  cm-2  at  750  °C,  1 .72  A  cm-2  at  700  °C,  and  1 .52  A  cm-2 
at  650  °C.  Fig.  6  shows  the  predicted  anode  concentration  polariza¬ 
tion  as  a  function  of  current  density  at  700  and  800  °C  calculated 
using  the  anode  limiting  current  densities  and  Eq.  (12).  The  anodic 
concentration  overpotentials  shown  in  the  figure  are  a  little  higher 
than  those  we  have  reported  previously,  which  is  due  to  the  fact  that 


In  this  work,  the  effective  binary  diffusivity  of  FI2-H2O  in  a 
porous  Ni-YSZ  anode  was  measured  using  an  electrochemical  cell 
under  the  operating  conditions  of  SOFCs.  The  effective  binary  dif¬ 
fusivity  of  FI2-FI2O  was  obtained  from  the  relationship  between 
the  applied  current  in  the  oxygen  pump  and  H2  partial  pressure 
gradient  across  the  porous  Ni-YSZ  anode  sample  measured  using 
the  oxygen  sensor.  The  FI2-FI2O  effective  binary  diffusivity  values 
ranged  from  0.070  cm2  s_1  at  800  °C  0.060  cm2  s_1  to  0.046  cm2  s_1 
at  650  °C.  The  anode  limiting  current  density  and  concentration 
polarization  were  estimated  using  the  experimental  results.  Future 
work  will  involve  improving  the  anode  microstructure  to  minimize 
the  concentration  polarization  using  the  technique  presented  in 
this  paper. 
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